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We examine the internal structure of solar-like stars in detail between 0.8 and 1.4 Mq and 
during pre-main sequence phase. Recent opacity computations of OPAL along with a new hy- 
drodynamical mixing process have been considered. We also introduce up-to-date nuclear reac- 
tion rates and explore the impact of accretion, mixing-length parameter, non-solar distributions 
among metals and realistic rotation history. We compare models predictions of lithium deple- 
tion to the ''Li content observations of the Sun and to 4 young clusters of different metallicities 
and ages. We show that we can distinguish two phases in lithium depletion: 1- a rapid nuclear 
destruction in the T-Tauri phase before 20 Myrs : this is independent of the mass used within 
our range but largely dependent on the extension and temperature of the convective zone, 2- a 
second phase where the destruction is slow and moderate and which is largely dependent on the 
(magneto)hydrodynamic instability located at the base of the convective zone. 

In terms of composition, we show the interest on considering helium and especially the mixture 
of heavy elements : carbon, oxygen, silicium and iron. We outline the importance of 0/Fe ratio. 
We note a reasonable agreement on lithium depletion for the two best known cases, the Sun 
and the Hyades cluster for solar-like stars. Other clusters suggest that processes which may 
partly inhibit the predicted premainsequence depletion cannot be excluded, in particular for 
stars below ~ O.QM©. Finally we suggest different research areas such as initial stellar models 
and more realistic atmospheres which could contribute to a better understanding of this early 
phase of evolution and which should become the object of subsequent research. 

Subject headings: Stars: pre-main sequence, abundances, interiors, rotation — Convection 



ABSTRACT 



Dynamical effects have been mainly ignored in 
classical stellar evolution during several decades 
even if they have been explored theoretically 
(Zahn 1974, 1992). Nowadays, helioseismology 
provides observational constrains on such effects 
and therefore allow us to begin to introduce them 
in the description of main sequence stars (Frolich 
et al. 1997; Gabriel et al. 1997; Kosovichev et al. 
1997). It permits observations of the solar convec- 
tive layers complex motions. Moreover acoustic 
mode determination allows the extraction of the 
rotation and sound speed profiles down to the 
energy-generation core (Kosovichev et al. 1997; 
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INTRODUCTION 



Dziembowski 1998; Turck-Chieze et al. 1997). 
Meridional circulation begins to be accessible to 
the solar seismic observations and some (magneto) 
hydrodynamical instability has been put in evi- 
dence at the base of the convective zone. In order 
to see the consequences of such a process, we have 
focused our attention on two elements ''Li and 
^Be which are destroyed in stars from the center 
to regions located slightly below the natural tran- 
sition between transport of energy by radiation 
and by convection. Lithium surface abundances 
history is directly related to physical processes 
in this region which is at the same time impor- 
tant for understanding dynamos in stars and also 
probably the role of the internal magnetic field in 



general. Brun, Turck-Chicze, & Zahn (1999) have 
invoked the hydrodynamical instability proposed 
by Spielgel & Zahn (1992) to explain both recent 
hclioscismic results and the solar lithium deple- 
tion. They introduced a turbulent term in the 
equation of diffusion. This term is directly related 
to local rotation and differential rotation. Evolu- 
tion of such a mixing term with recorded surface 
rotation history of sun-like stars in open-clusters 
finally explains the present observed solar photo- 
spheric element abundances. It also produces the 
right order of magnitude of lithium destruction as 
a function of time for open-clusters older than a 
Gyr (NGC 752, M67). 

In this paper we focus on the different in- 
gredients of the mainsequencc (MS) and above 
all premainsequence (pre-MS) structural evolution 
of solar-like stars through the excellent indicator 
which is the ''Li sin-faco abundance. From that 
viewpoint stellar evolution models have received 
a lot of attention for a long time. Almost forty 
years ago pro-MS evolution was already proposed 
to explain the low solar Li abundances relative 
to the solar system value as well as field stars 
and Hyades depletion pattern with temperature 
(Bodenheimer 1965). Bodenheimer described the 
main features of lithium burning in pre-MS. Accu- 
rate ^Li observational data for open clusters have 
similarly been obtained for long (Zappala 1972). 
They have outlined the complex history of this el- 
ement and suggested depletion with age on MS. 
In the past twenty years many theoretical and ob- 
servational works have addressed the topic. These 
investigations suggested that the understanding of 
lithium abundance might be related to many dif- 
ferent non-standard processes extending from mi- 
croscopic diffusion (Michaud 1986) to large scale 
mixing through rotation (Baglin et al. 1985) an- 
gular momentum evolution (Pinsonneault et al. 
1989) or internal waves (Schatzmann 1993). Dur- 
ing pre-MS lithium depletion strongly depends on 
typical temperatures within convection zones and 
therefore input physics (Profitt & Michaud 1989). 
Until today there is no safe explanation of lithium 
history among solar-like stars. It is very likely that 
various phenomena are indeed responsible for pre- 
MS and MS Li evolution as well as G-type stars 
depletion pattern and lithium dip. This contro- 
versial subject therefore remains a broad and ac- 
tive area of research (D'Antona & Mazzitelli 1994, 



1997; Ventura et al., 1998; Palla 1999; D'Antona, 
Ventura & Mazzitelh, 2000). 

This work intends to show results involving re- 
cent opacity and nuclear reaction rates determina- 
tion. We consider possible variations of composi- 
tion and non solar repartition among metals as in 
the next decade many improvements in the knowl- 
edge of the detailed photospheric composition of 
the young stars are anticipated. We envision vari- 
ations of mixing-length parameter with age and a 
new hydrodynamical instability along with a real- 
istic rotation evolution. We outline the impact of 
such improvements on the topic of early phases of 
evolution which is meant to be better constrained 
by future seismic informations. In section 2, we 
describe the bases of the stellar models we use and 
examine the solar case in details. The role of the 
different elements through their composition and 
the opacity coefficients is discussed in details in 
section 3. Section 4 is dedicated to convection & 
accretion and section 5 to rotation. Finally we 
give some perspectives in section 6. 

2. PRE-MAIN SEQUENCE EVOLUTION 
OF SOLAR-LIKE STARS 

2.1. Physics and global evolution 

In this study, we consider models in the mass 

range from 0.8 to 1.4 Mq and different composi- 
tions corresponding to the Sun, the Pleiades, and 
the Hyades. Our starting point is a fully adia- 
batic polytrope with a central temperature of ~ 
3 10^ K and a radius of 20 Rq for IM© star. So 
the computation begins prior to deuterium burn- 
ing which corresponds to the observed birth-line 
(Stabler 1988). We observe that models are still 
fully convective polytropes when the star reaches 
the birth-line with a radius of about 6 Rq so we 
expect no consequences of this early-phase com- 
putation except the idea to take as an hypothesis 
an initial polytropic structure. 

Models have been computed using the CESAM 
code (Morel 1997) and all the updated physics use- 
ful for the refined solar models (see Brun, Turck- 
Chieze, & Morel 1997). For the reaction rates of 
the pp chain and CNO cycle, we consider the com- 
pilation of Adelberger et al. (1998) and for the 
''Li{p,a)^He the work of Engstler et al. (1992). 
We use OPAL equation of state (Rogers, Swen- 
son, & Iglesias 1996) and opacities (Iglesias & 
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Rogers 1996) above 5800K. For lower tempera- 
tures the OPAL equation of state is replaced by 
the MHD equation of state (Mihalas, Dappen, & 
Hummer 1988) and the opacities arc from Alexan- 
der & Ferguson (1994). Interpolation of opaci- 
ties are performed using v9 birational spline pack- 
age of Houdck (Houdck & Rogl 1996). The at- 
mosphere is connected to the envelope at optical 
depth 10 where the diffusion approximation for ra- 
diative transfer becomes valid (Morel ct al. 1994). 
The opacity atmospheric is Rosseland mean opac- 
ity extracted from Alexander & Ferguson (1994) or 
Kurucz (1992) sets. We therefore do not use grey 
approximation. Surface boundary conditions are 
p = 3.55 10~^g.cm~^ when r = 10~^. Standard 
mixing-length theory is applied and the convec- 
tion zone is completely chemically homogeneous. 

During the pre-MS, stars are in quasi- hydrostatic 
equilibrium and slowly contract towards ZAMS 
within a time scale comparable to their Kelvin- 
Helmoltz time-scale. Defining ZAMS as the age 
where the thermonuclear hydrogen fusion provides 
99% of stellar energy the pre-MS lifetime varies 
from 30 Myrs for a star of 1.4 Mq Pleiades com- 
position up to 100 Myrs for a 0.8 Mq star of the 
Hyades composition. The Sun lies in between, 
with a pre-MS of ^ 50 Myrs. The transition from 
fully convective object to radiative core structure 
depends on mass and composition. Higher mass 
stars contract faster, increase internal temper- 
atures faster and so decrease radiative thermal 
gradients faster. Moreover lower metallicity accel- 
erates contraction and decreases opacities which 
also favor radiative stratification. Therefore lower 
metallicities and higher masses give rapid raise to 
radiative core. A massive (1.4 Mq) Pleiades com- 
position star develops radiative core after 0.85 
Myrs and a 0.8 Mq Hyades composition star after 
3.5 Myrs. In the solar case, radiative core exceeds 
1% of total mass at 1.8 Myrs and then fastly ac- 
celerates in mass and radius to reach dimensions 
very close to its actual dimensions at 25 Myrs. 

Due to the swift convective movements, surface 
matter is repeatedly exposed to physical condi- 
tions that prevail in deep interior. Figure 1 illus- 
trates this point in showing the time dependent 
evolution of the thermodynamical quantities (tem- 
perature and density) at the base of the convection 
zone (hereafter BCZ) in regard to the photospheric 
lithium for the case of a young Sun. At the be- 
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Fig. 1. — Time evolution of the temperature (con- 
tinuous line) and density (dashed line) at BCZ for 
a 1 Mq star of solar composition. Both tempera- 
ture and density are normalized to their maximum 
values, respectively 3.94 x 10® K and 1.97g/cm'^. 
We have also mentioned the density dependence of 
the temperature for a characteristic burning time 
of 1 Myrs (long-dashed line) and the lithium pho- 
tospheric fraction (dot-dashed line) normalised to 
maximal value 9.3 x 10~® in mass fraction (ie 3.2 
dex). 
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ginning of the considered evolution BCZ (which 
coincides with the stellar center) is too cool to al- 
low lithium burning. As the star evolves on pre- 
MS, deep regions of convection zone temporally 
exceeds "^Li burning point in typical stellar condi- 
tions (~ 2.510^K). "^Li therefore offers a direct 
insight over stellar internal structure and evolu- 
tion as it is extremely sensitive to the appearance 
of the radiative core. We note that early "^Li de- 
pletion occurs in the very beginning of pre-MS. Ef- 
fectively, for IMq, the BCZ temperature increases 
rapidly from less than 10^ K to approximately 4 
10® K at 2 Myrs, the density increases also (but 
not quite simultaneously) up to 2 g/cw? between 
7 and 8 Myrs. Then they slowly sink towards val- 
ues which are very near from the conditions of the 
BCZ of the present Sun. This evolution results in 
"^Li burning from 2 to 20 Myrs. The way de- 
pletion occurs presents several difficulties. First, 
depletion takes place just at the BCZ as convec- 
tion rapidly recedes in mass-fraction passing from 
the whole star mass at 1 Myr to less than 20% 
20 Myrs later. Moreover depletion typical time 
evolves rapidly and decays to very low values com- 
pared to time scale evolution of the stellar struc- 
ture down to ~ 10^ years at BCZ. These peculari- 
ties impose small time step and very thin meshes. 
Computations suggest that one third difficulty 
must arise. Early radiation zone stratification 
hardly differs from convection stratification. In 
other words the newly stabilized medium is (and 
stays for at least 10 Myrs) close to convective in- 
stability. We can illustrate this using polytropes : 
the adiabatic polytropc [Pap^'^^^'^ with n=3/2) is 
'harder' than its radiative counterpart (Pap^+^/" 
with n=3). Radiation zone which represents half 
of the stellar mass at 10 Myrs, should there be 
more concentrated in the core. Figure 2 shows the 
contrary. The fact that radiation energy trans- 
port replaces convective motions does not mean 
that radiative stratification immediately estab- 
lishes. At 10 or 20 Myrs the stratification in radi- 
ation zone is still very close to adiabatic stratifica- 
tion. This has consequences on rotation evolution 
as we shall see in section 5. We believe this prop- 
erty to be related to Kelvin-Helmoltz time value. 
The stellar interior needs at least this long to redis- 
tribute thermal energy and evolve from convective 
to radiative structure. This suggests that little 
modifications in stellar structure or new phenom- 



ena taken into account could easily change BCZ 
position at times from 10 to 20 Myrs. We may 
for instance suppose that even a low energy den- 
sity magnetic field would stabilize or destabilize 
the region. This point leaves a priori room for 
many non-standard mechanisms that will have to 
be investigated and that will probably lead to sub- 
stantial variations in pre-MS surface Li depletion 
(Ventura et al. 1998). 

2.2. Adjustment of evolution parameters 
for "^Li burning 

"^Li reaction rate is extremely sensitive to tem- 
perature. For typical prc-main sequence condi- 
tions, at the maximum temperature of BCZ (3.5 to 
4 10*^ K) this rate varies as T^^ to T'^^. The deple- 
tion is therefore very sensitive to conditions near 
the base of the convection zone (figure 3). Con- 
sequently, the photospheric Li depletion strongly 
depends not only on the mass of the convection 
zone but also on the precision of the mass shell di- 
vision at BCZ. This point is essential to perform 
valid integrations. We have adjusted the mesh in 
order to limit the variation of the reaction rate 
at BCZ between two layers to 15 %, this corre- 
sponds typically to meshes of 0.005 Mq or 0.003 
Rq. Moreover, the typical Li pre-MS burning 
requires also a very accurate resolution relative to 
time integration. Figure 3 shows the evolution of 
the characteristic time along the evolution. At 3- 
4 Myrs we find a characteristic destruction time 
at BCZ of only S.IO** to 1.310''' years (depending 
on composition). So we have adjusted the typical 
time step to be a factor 10 lower than the charac- 
teristic depletion time. In figure 3 we have shown 
three conditions corresponding to the maximum of 
the temperature at BCZ (3 Myrs), to the arrival 
on ZAMS and to the present age of the Sun. For 
each condition, the depletion time varies by ap- 
proximately one order of magnitude from BCZ to 
the right end of the plot but this extension corre- 
sponds to only 10 % variation of the mass of the 
convection zone. 

2.3. Solar lithium burning during pre- 
mainsequence 

Pre-MS Li depletion within stars of roughly 
solar mass and composition has been estimated 
in various previous investigations (Bodenheimer 
1965; D'Antona & Mazzitelh 1984 ; Profitt & 
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Fig. 2. — Solid lines present the mass distribution 
of polytrope vs radius for a IMq star. Upper solid 
line corresponds to a polytrope of index n=3 (rep- 
resentative of radiative stratification), lower solid 
line corresponds to a polytrope of index ti—3/2 
(adiabatic stratification). Dashed-lines present 
model stratifications at 10 Myrs (lower line) and 
30 Myrs (upper line) for solar composition and 
mass. Crosses on these lines are BCZ. It is note 
worthy that at 10 Myrs stratification is still almost 
adiabatic albeit more than 50% of stellar mass is 
in the radiation zone. 
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Fig. 3. Temperature dependence of the Li de- 
pletion characteristic time around BCZ for a dif- 
fusive solar model (X=0.70821, Y=0.2722). The 
continuous line corresponds to the solar age, 
dashed line to the ZAMS position, and dot-dashed 
line to a 3 Myrs star. Triangle, cross, and diamond 
represent the positions of the BCZs. 



Table 1: Ratio of lithium fraction to initial frac- 
tion after pre-MS depletion phase within solar 
composition stars (see text for details). Our re- 
sults (PTCOl) are compared to previous ones of 
D'antona & MazziteUi 1984 (D'A&M84), Profitt & 
Michaud 1989 (P&M89) and Ventura et al. 1998 
(V98). P&M89a and P&M89b respectively cor- 
responds to Z=0.0169 and 0.024. For l.l&lM© 
no depletion is predicted on MS. For O.9M0 star 
D'antona & Mazzitelli do not predict depletion on 
MS but Profitt & Michaud 1989 and us predict it. 
For these low mass stars lithium fraction therefore 
refers to a fixed age of 70 Myrs. Results of Ventura 
et al. 1998 are provided here only for solar mass 
stars as they do not give results for other masses 
in MLT framework. 



Mass (Mq) 


0.9 


1 


1.1 


D'A&M84 


0.664 


0.871 


0.949 


P&M89a 


0.467 


0.679 


0.823 


P&M89b 


0.118 


0.317 


0.535 


V98 




0.026 




PTCOl 


0.002 


0.071 


0.244 



Michaud 1989; Ventura et al. 1998). These three 
last studies suggest lithium depletion to increase 
as input physics is updated. D'Antona & Mazz- 
itelli 1984 find that a solar mass star depletes 
~ 13% of its initial lithium. The authors choose 
Y=0.23, Z=0.02 and amit = 2. Profitt & Michaud 
1989 find that a solar mass star depletes from 
~ 32% to ~ 69% of its initial lithium. These 
values correspond to the same amit = 1-5 but 
make different assumptions on the composition. 
The first one refers to Y=0.28, Z=0.0169 as the 
second one refers to Y=0.28, Z=0.024. Profitt 
& Michaud suggest that difference to D'Antona 
& Mazzitelli depletion rates is probably due to 
opacities improvements. They moreover remark 
that htc gap would get even wider if they were 
to use D'Antona & Mazzitelli parameters (Y and 
amit)- The gap between 1984 and 1989 results in- 
creases towards low masses and gets even wider 
if one considers very recent work from Ventura et 
al. (1998). A model from these authors having 
solar mass and composition (Y=0.28 and Z=0.02) 
and using MLT amit = 1-55 brings initial lithium 
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Fig. 4. Time dependence of the photospheric 
Li abundance for a IMq star with solar initial 
composition for different computations: a diffu- 
sive calibrated model (Y=0.2722 and Z=0.01959) 
with a time step adjusted to Li BCZ burning time 
(continuous line), idem but the time step imposed 
by the stmcture of the star (dot-dashed line), a 
nondiffusive noncalibrated model with a time step 
adjusted to "^Li BCZ burning time (dashed-line), 
idem but the time step imposed by the structure; of 
the star (long-dashed line). The calibrated nondif- 
fusive model with Y=0.2648 and Z=0.01763 with 
correct time step is represented by (-...-...) line 
and shows the effect of the composition 

fraction 3.3 dex down to 1.72 dex (a decrease of 
a factor ~ 40). Altough a bit larger, such deple- 
tion is the same order of magnitude as what our 
computations suggest. The present work therefore 
confirms the tendency of depletion to increase with 
updated physics. We find that lithium lessens by 
roughly a factor ten in solar mass and composition 
star (we take Y=0.2722, Z=0.01959, a™;, = 1.766). 
Opacities must partly be responsible for this evolu- 
tion but new estimations of nuclear reaction rates 
should also be responsible for this increase. The 
Engstler et al. (1992) Li{p,aYHe astrophysical 
factor is indeed increased by 30 % in comparison 
with older values (see BTZ99). Finally the correct 
attention to the adapted timestep to lithium de- 
pletion rate is a third possible source of dispersion 
(as shown below). Table 1 compares our compu- 
tations to previous ones. 

Figure 4 presents the results for the evolution of 



the Sun in pre-mainsequence. Following the work 
of Brun, Turck-Chieze, & Zahn (1999, hereafter 
BTCZ99) we have computed models of the Sun, in 
introducing (or not introducing) microscopic diffu- 
sion. Depending on this choice, the initial compo- 
sition Y= 0.2722, Z= 0.01959 (or Y= 0.2624, Z= 
0.01763) is adjusted to get the correct luminos- 
ity and radius at the present age, calibrated with 
an accuracy better than 4 10~^. Figure 4 shows 
that the microscopic diffusion has no effect on the 
pre-mainsequence. It is too slow a process as it re- 
sults in a change of photospheric helium and met- 
als abundances of only 10% along the whole life of 
the Sun (Turcotte et al. 1998). We have to note 
here a crucial point in ''Li pre-MS depletion. Sun- 
like stars change much faster on pre-MS than on 
MS because evolutionary time scale is primarily 
related to contraction (hence to Kelvin-Helmoltz 
time tkh) which is much smaller than nuclear 
reaction time (tat). Yet one should not expect 
any calculation on lithium during early pre-MS 
to provide correct results without some cautions. 
Despite tkh « tn forces to adopt in the code 
small temporal steps these are much larger than 
^Li burning time at BCZ (rrj^j) at these ages. The 
following comparisons illustrates this. In diffusive 
models ''Li depletion is increased by a factor of 
3.5 between a star where time evolution is always 
less than one tenth of and a model where the 
time step is given by the evolution of the struc- 
ture. For nondiffusive models this ratio goes up to 
5.3. Once time step is well below rr^^j no change 
in ''Li depletion is seen at a given composition for 
various timesteps and with or without diffusion. 
In every following computations we take time step 
to be as long as lithium is rapidily depleted 
(prior to ~ 30 Myrs). Furthermore burning rate 
is mass-averaged over convection zone. We have 
checked that the results are robust to the choice of 
time step. If we decrease time step down to ^^^^ : 
''Li depletion increases of only 3.5 %. With time 
step adjusted for a good treatment of 'Li burn- 
ing, the difference of about 0.5 dex after 20 Myrs 
(a factor 3 in destruction) between diffusive and 
nondiffusive models is only due to the change of 
initial composition. This shows the great depen- 
dence of the lithium burning on the composition 
(see next section). Such precise computation of 
the lithium burning in the pre-mainsequence was 
not included in the previous work of BTCZ99. 
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Of course the notion of calibration between dif- 
fusive and nondiffusive models is not justified in 
the study of young clusters but the great sensitiv- 
ity of this calibration shows already the difficulty 
of predicting lithium burning in pre-mainsequence 
stage. 

Li is generally not the only probe of the stel- 
lar internal structure but ^Li is depleted at lower 
temperature and is of no help here. We have ob- 
served that Be burning gives similarly poor indi- 
cations. In fact, computed stellar models of 0.8 to 
1.4 Mq show a ^Be depletion of less than 0.04 dex 
whatever the cluster membership (from Pleiades 
to Hyades). 

3. THE ROLE OF THE DETAILED COM- 
POSITION 

3.1. Open cluster evolution 

Open clusters allow a direct test of ''Li evolu- 
tion with time and/or composition. In this pa- 
per we focus our attention on two young and 
two midle-age clusters : Pleiades and Blanco I 
{C, Sculptoris) for the former Hyades and Coma 
Berenices for the latter. The objective is to dis- 
tinguish mctallicity from age effects. Blanco I and 
the Hyades apparently are metal-rich clusters with 
respectively [Fe/H] = 0.127±0.022 (Boesgaard & 
Friel 1990) and [Fe/H] = 0.14 ± 0.01 (Jeffries & 
James 1999) however Blanco I is much younger 
than the Hyades and with an age around 50 to 
90 Myrs comparable to the Pleiades cluster which 
is estimated to be 70 Myrs old (Patenaude et al. 
1978). On the other hand Pleiades and Coma 
Berenices clusters have mctallicities very slightly 
below solar one with [Fe/H] = -0.034 ± 0.024 
(Boesgaard & Priel 1990) and [Fe/H] = -0.052 ± 
0.026 (Friel & Boesgaard 1992) respectively but 
Coma has an age of 500 Myrs comparable in age 
to the old Hyades cluster of ~ GOOMyrs (Perry- 
man et al. 1998). It should be mentioned that 
composition data are undoubtedly more reliable 
in the case of Hyades or Pleiades than in the case 
of Blanco I which is more distant (^ 250 ± 30pc) 
than all other clusters and has been less studied. 

Unlike the sun, microscopic diffusion has a neg- 
ligible effect in the evolution of these young or 
midle-age clusters. We compute a decrease of 1.3 
and 0.8 % for helium and metals in Hyades solar 
mass star case. As such variations always remain 



largely smaller than error bars over mctallicity or 
helium content it is immaterial to study initial 
composition effect resulting from diffusion. This 
is particularly true for lithium. Let us however 
remark that microscopic diffusion time appears 
roughly ten times shorter for F type stars. We 
compute a decrease of respectively 10% and 6% 
for helium and metals in 1.4Mq {T^ff ~ 66001^) 
at Hyades age. However at an age of 50 Myrs 
decrease of helium and metals arc only 0.6% and 
0.4% so that it is not plausible that microscopic 
diffusion changes early lithium history even in case 
of slightly more massive stars than the sun. 

In this study, we do not consider different ini- 
tial "^Li abundances. In all four clusters stars 
more massive than 1.4 solar masses (Teff>6900 
K) exhibit the same ^ii content of 3.2 to 3.3 dex 
for Pleiades, Hyades, (Soderblom et al. 1993) 
and Coma (Boesgaard 1987) . Blanco I depletion 
pattern is furthermore identical to Pleiades' one 
(Jeffries 1999). The early F-stars abundances re- 
main unchanged and arc compatible both with 
very young T Tauri stars (Magazzu, Rebolo, 
& Pavlenko 1992) and the interstellar medium 
present '^Li value (Knauth et al. 2000). Indeed, 
it seems that there was no significant evolution in 
galactic gas ^Lz/hydrogen ratio over the last 1.7 
Gyr (Hobbs ct al. 1988) so that the initial ''Li 
does not vary from Hyades formation time until 
today. For all the clusters we take the same stan- 
dard value of 3.27 dex for initial ''Li abundances. 

In the following we will investigate the effects of 
the composition on the lithium burning in separat- 
ing the effects of deuterium, helium, and metals. 



3.2. Sensitivity of the lithium burning to 
the deuterium composition 

Deuterium, the most fragile element, is de- 
pleted around 5 10^ K in stellar interiors. The 
galactical evolution leads to a continuously de- 
crease because of astration. Recent measurements 
show an abundance of {D/H)ism = 1.46 ±0.09 x 
10-5 (Piskunov et al. 1997) or {D/H)ism = 
1.60 ± 0.09 X 10-5 in direction of CapcUa (Linsky 
et al. 1995) whereas presolar value is estimated 
as {D/H)preQ = 3.01 ± 0.17 X 10-^ by Gautier & 
Morel (1997). Young clusters initial abundances 
are probably the same as present ISM one. 
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During contraction, the stars experience sev- 
eral stages of light element burning. First, deu- 
terium burning (T = 5 10^ K) significantly partic- 
ipates to the energy production. Palla & Stabler 
(1991) have shown that for 1 Mq star, deuterium 
burning stops the contraction at a radius of ~ 
5 — 6i?0 which defines the so called birth-line 
(Stabler 1988). The deuterium burning stops be- 
fore 1 Myr when first phase lithium depletion 
starts. Consequently the star should have enough 
time to 'forget' its previous history. We checked 
this briefly in the framework of our models. For 
solar mass stars, deuterium begins to burn at the 
age ~ 4 10^ years during ~ 2 10^ years duration. 
Then lithium depletion starts after 1.4 10^ years 
shortly prior to radiative core appearence. We 
have considered three deuterium mass fraction: 
the solar value, the actual ISM value and also a 
null value. Such variations hardly change lithium 
photospheric abundances evolution. Models hav- 
ing Hyades or Pleiades composition, show maxi- 
mum variations between them of ~ 10% (0.04 dex) 
up to clusters or solar ages. When deuterium frac- 
tion decreases, lithium depletion slightly increases. 

3.3. Sensitivity of lithium burning to the 
heUum content 

"^Li evolution is sensitive to metals but also to 
helium mass fraction. The solar photospheric he- 
lium has been determined by helioseismology to be 
0.246-0.249 (Basu & Antia 95). Including micro- 
scopic diffusion (partly inhibited by turbulence) in 
calibrated solar models provides a mean to reach 
the initial solar helium content, estimated to be 
0.2722 (BTCZ99). For open clusters there is no 
recent helium abundance determinations for any 
studied cluster except the Hyades. A ZAMS low 
mass stars fit of Hyades suggests the cluster helium 
fraction to be 0.26 ± 0.02 (Ferryman et al. 98). 
This value is confirmed by the position of a num- 
ber of Hyades binaries in M-L diagram (Lebreton 
2000). We adopt this value in our calculations. 
However other calculations (Pinsonneault et al. 
1998) suggest an higher value of 0.283. For the 
other clusters we use the usual helium to metallic- 
ity scaling law : AY/AZ=3 ± 2 as deduced from 
calibration of nearby visual binary systems (Fer- 
nandes et al. 1998). In the same manner eval- 
uations based on HII regions suggest an increase 
of helium content with metals. If we refer to fig- 



ure 12 from Meyer 1989 we find that AY/AZ lies 
between 1.5 and 6. Such a trend is not totally 
satisfactory as the Hyades do not seem to follow 
this law. In fact the use of such a law requires 
a knowledge of both helium and metallicity in at 
least one site which is generally taken to be our 
solar system neighborhood. But the solar system 
seems not to be representative of 4.6 Gyrs ago 
mean trends of the ISM (Gies & Lambert 1992) 
so it is possible that deduced helium fractions are 
biased. Other calibrations of the metal/helium re- 
lation coming from globular clusters and galactic 
bulge measurements do not lead to similar helium 
fractions (Deliyannis, Demarque & Kawaler 1990). 
Changes here are not so important, so we retain 
here solar systcmi abundances as reference point. 

We have estimated the impact of helium vari- 
ation on lithium depletion around solar compo- 
sition, for Pleiades age, together with reasonable 
variations of the mixing-length parameter a, and 
metal fraction. The effect of the mixing length pa- 
rameter is small, an increase of 5% leads to a de- 
crease of "^Li content by 20%. On the contrary the 
impact of the composition is high: an increase of 
A Y=: 0.025, leads to a decrease of lithium burning 
by 64%, the corresponding A Z = 0.025/3 leads to 
an increase of lithium burning by a factor 300 and 
the correlated variation of helium and metallicity 
leads to an increase of lithium burning by a factor 
60. We address the metallicity dependence more 
precisely in the next chapter. The anti-correlation 
between lithium depletion and helium is reported 
in previous calculations. D'Antona & Mazzitelli 
(1984) analyse hefium mass fraction 0.23 and 0.28. 
There the rate of lithium fraction remaining from 
their initial fraction to this helium mass variation 
i^/AY) is roughly 3.3 and 1.3 for 0.9 and IMq 
stars respectively. Considering a helium variation 
(~ 0.05) from 0.2624 to 0.32 we correspondingly 
find rates 2.7 and 4.7 for 0.9 and IMq stars. This 
is the same order of magnitude as D'Antona & 
Mazzitelli one although our dependence seems a 
bit larger. The sensitiveness on helium fraction 
can easily be understood since opacities in stellar 
interiors reduce with helium content : less elec- 
trons are available for the same amount of mat- 
ter. If we assume that helium content is a free 
parameter we can speculate as to what level it 
should be increased to agree with the observations 
in Pleiades case. Lithium abundances being scat- 
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tered for any given effective temperature one has 
first to derive a mean value. Figure 5 (daslied line) 
shows a third degree polynomial least square fit of 
observed abundances. Wc obtain helium fraction 
of 0.36 for 0.85Mq star, 0.32 for IMq, and 0.3 for 
IAMq. These very high values seem difficult to 
justify. Recent work by Dcharveng et al.(2000) on 
helium abundances in HII galactic regions shows 
He/H to be below 0.105 where it can safely be de- 
termined from He^/H^ ratios. This corresponds 
to helium mass fraction below 0.3. Another work 
indicates helium mean mass fraction should be 
0.28 ± 0.02 in the galactic bulge (Minniti 1995) 
which is probably an upper limit to galactic abun- 
dances as most chemically evolved stars popula- 
tions are expected to be near the center of our 
galaxy. A second difficulty is that helium mass 
fraction would have to vary with mass by about 
20%. This certainly is not what we observe in 
Pleiades. Nevertheless, it is clear from this analy- 
sis that a proper determination of helium is a cru- 
cial knowledge to understand the young clusters 
lithium evolution. 

3.4. Sensitivity of lithium burning to the 

metallicity 

We have computed several models to describe 
the 1 Mq Hyades and Pleiades stars. These mod- 
els rely on physics described in section 2.1. Ta- 
ble 2 summarizes these results. We have calcu- 
lated three types of models: the first one (A) uses 
metallicity deduced from the observation of iron 
given in section 3.1 and a parameter calibrated 
on solar models (1.766). The second model (B) 
adopts the same composition but three different 
values of a were used to account for recent 2D hy- 
drodynamical evaluations (see section 4.1). The 
third model (C) is an extremal model in the sense 
that it adopts the lowest metallicity and highest 
helium content within error bars in an attempt 
to cancel discrepancies between calculations and 
observations. In the Hyades case we have added 
a supplementary model (D) introducing the value 
of 0.283 in helium mass fraction claimed by Pin- 
sonneault et al. 1998. Other Hyades models have 
been computed, we discuss them in more detail in 
section 3.5 

On the main sequence effective stellar temper- 
atures evolve slowly. The temperature of an open 
cluster solar mass star will therefore be reliable in- 



dependently on age uncertainties. Moreover these 
temperatures do not vary much within composi- 
tion uncertainties or a parameter different evalu- 
ation along early pre-main sequence (sec table 2). 
For instance the IMq Hyades models exhibits a 
variation of less than 10 K around 5450 K between 
550 and 700 Myrs (idem for the model including 
a parameter effects on early pre-MS variation but 
around 5500 K). 

We note in figures 5 and 6 that ^Li depletion 
is too strong on pre-main sequence for solar-type 
stars with Pleiades or Hyades composition. This 
agrees with recent results from Morel et al. (2000) 
of the B component of l Pegasi binary system. 
With an estimated age of 56 Myrs for the sys- 
tem and [Li]=2.69 dcx for the B component this 
0.819 star is clearly underdepleted. This problem 
is pointed out in other recent studies. Ventura et 
al. (1998) using the full spectrum of turbulence of 
Canuto, Goldman, & Mazzitelli (1996) prescrip- 
tion for modelling the convection found also a too 
strong depletion for solar composition stars, and 
similarly a very strong dependence on metallicity 
and mass. 

It is well known that open cluster stars exihibit 
an anti-correlation between effective temperature 
and lithium abundances. Moreover the dispersion 
in lithium abundances grows when temperature 
declines. This dispersion is too large to be due to 
abundance uncertainties (Soderblom et al. 1993). 
This is also too large to be due to color errors as 
the vector of temperature errors is nearly paral- 
lel to the mean lithium- Tg/j trend in cool open 
cluster stars. Moreover this dispersion is already 
observed in ^Li equivalent width as a function of 
color (Soderblom et al. 1993, Thorburn et al. 
1993). At a given effective temperature there is 
undoudtedly real star to star differences. 

Observations for each star give an effective tem- 
perature and a lithium abundance but these quan- 
tities are not directly measured as they are de- 
duced from photometric measurements such as (B- 
V) and equivalent width of a particular spectral 
absorption line. Measurements are always slightly 
scattered and moreover there are unavoidable er- 
rors when doing conversions. For Pleiades open- 
cluster, Soderblom et al. (1993) find uncertainty 
for effective temperature of ~ 130K and uncer- 
tainty for "^Li ~ 0.05 dex. For the Blanco I 
cluster Jeffries & James (1999) give an effective 
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Table 2: IMq stellar models for Sun, Pleiades, and Hyades compositions. All models include microscopic 
diffusion. Several a parameter values have been considered : 1.766 is value assumed at solar age it is deduced 
from solar calibration. a= 1.748 results from slightly different calibration when tachocline mixing is taken into 
account. 1.935 and 1.850 are values induced from hydrodynamical simulations (Ludwig, Preytag, feSteffen 
1999). For Pleiades composition we take Q!=1.935 before 15 Myrs and a=1.85 from 15 to 22 Myrs. For 
Hyades composition we take q;=1.935 before 20 Myrs and a=1.85 from 20 to 28 Myrs. Lithium dex fraction 
is presented at appropriate age. The last column provides effective temperature to outline composition 
impact on it. 





Z 


Y 




a 


^Li content 


[0/Fe] 


Tachocline 


'^eff 




(mass fraction) 


(mass fraction) 






(dex) 






(K) 


Solar 


















models 


















Reference 


1.959 10-^ 


0.2722 




1.766 


2.0 





No 


5776 K 


Tachocline 1 


1.903 10-^ 


0.2695 




1.748 


0.35 





Yes Tdisk = 0.5 


5777 K 


Tachocline 2 


1.903 10-^ 


0.2695 




1.748 


1.1 





Yes Tdisk = 10 


5777 K 


Pleiades 


















models 


















A 


1.632 10"^ 


0.2624 




1.766 


2.6 





No 


5717 K 


B 


1.632 10"^ 


0.2624 


1.935, 


1.850, 1.766 


2.4 





No 


5715 K 


C 


1.535 10"^ 


0.2679 




1.766 


2.75 





No 


5801 K 


Hyades 
models 


A 


2.367 10"^ 


0.2633 




1.766 


0.9 





No 


5450 K 


B 


2.367 10-^ 


0.2633 


1.935, 


1.850, 1.766 


0.5 





No 


5466 K 


C 


2.180 10"^ 


0.28 




1.766 


1.75 





No 


5646 K 


D 


2.28 IQ-^ 


0.283 




1.766 


1.55 





No 


5632 K 


E 


1.57 10-^ 


0.26 




1.766 


2.64 


-0.2 


No 


5802 K 


F 


1.94 10"^ 


0.26 




1.766 


1.98 


-0.2 


No 


5638 K 


G 


1.57 IQ-^ 


0.26 




1.766 


2.3 


-0.2 


Yes Tdisk = 10 


5805 K 


H 


1.57 10-^ 


0.26 




1.766 


1.75 


-0.2 


Yes Tdisk = 0.5 


5804 K 



temperature uncertainty of ~ 250K and an abun- 
dance uncertainty of 0.11 dex. Figures 5 and 6 
clearly show that such uncertainties are not large 
enough to recover agreement to computations. As 
we mainly analyze IMQ evolution, we concen- 
trate on corresponding error-bars. It is gener- 
ally found in litterature that these stars verify 
0.6 < (B - V) < 0.75. In fact, lithium abun- 
dances are such steep functions of effective tem- 
perature that if we use this criterium, we gather 
in the same sample stars with very different "^Li 
fractions. Let us consider Hyades case. In Thor- 
burn et al. (1993) data, the mean value of lithium 
content for stars exihibiting 0.6 < {B — V) < 0.75 



is NCLi) = 2.0 ± 0.6 dex, but ii we limit to 

5400if < Teff < 5500K which is what we get 
in computation, we get N{J Li) = 1.3 ± 0.1 in the 
same data set. The small sample presently dis- 
cussed greatly reduces differences between obser- 
vations and numerical predictions. Yet differences 
are not definitively cancelled out, and it is obvious 
from figure 5 and 6 that such possible misleading 
associations between stellar mass and (B-V) will 
not succeed in explaining the very large observa- 
tion/theory gap concerning lower mass stars. 
Table 2 shows that expected dispersions among 
metallicity or uncertainties on helium fraction can 
be responsible for dispersion of lithium abun- 
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dances. In Pleiades case, the metallicity variation 
from -0.034 dex (Z= 1.632 lO'^) to -0.058 dex 
(Z= 1.535 10~^) minimal value corresponds to 6% 
in metal mass fraction and produces 0.35 dex vari- 
ation in photospheric lithium, which is quite the 
order of observed dispersion for solar mass stars 
within this cluster. In Hyades case, helium uncer- 
tainties at a level of ~ 7% correspondingly pro- 
duce 0.6 dex "^Li dispersion between models types 
(A) and (D). Composition variations are therefore 
able to explain ^Li dispersion through pre-MS ini- 
tial depletion effect but mass accretion could con- 
tribute to it also (section 4). If such variations 
explain lithium dispersions in young Pleiades or 
Coma clusters they also have to deal with middle- 
aged cluster observations. In Hyades dispersion in 
the Tf.ff—'^ Li relation has considerably decreased. 

3.5. Opacity role on "^Li burning 

Because of strong dependence of lithium pre- 
MS depletion on metallicity opacities effects 
in stellar interiors have been studied (Turck- 
Chieze et al., 1993; Turck-Chiczc 1998; Turcotte 
& Christensen-Dalsgaard 1998). Solar- like star 
metal opacities increase is responsible for the tran- 
sition between radiative and convectivc energy 
transport. In the present Sun, the main contrib- 
utors to the opacity at BCZ are oxygen and iron. 
They correspond respectively to 36 and 20 % of 
the total opacity (table 3). Situation is somewhat 
similar regarding pre-main sequence but the prin- 
cipal metal contributors to opacity are not the 
same because the medium is denser and hotter. 
The roles of oxygen and iron are of the same order 
of 20 %, then neon, silicium, and magnesium each 
represent approximately 10 % of total opacity. 
Up to now we have deduced metallicity from 
[Fe/H], in keeping the solar distribution inside the 
metals. There have been hints that solar abun- 
dances may deviate from mean ISM abundances 
and that solar system would not be representative 
of general trends (Gies & Lambert 1991). Regard- 
ing open clusters, similar situation could appear 
even if direct determinations have to be consid- 
ered cautiously due to chromospheric activity see 
Cayrel et al. (1985) for Hyades. Let us consider 
Hyades cluster. It is the closest and probably 
the best documented. Moreover it appears metal- 
rich. Hence Hyades stars experience amplified 
early pre-MS "^Li depletion. As we have already 
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Fig. 5. — Crosses : Li abundances for members 
of the Pleiades (Soderblom 93 data). Diamonds : 
Coma Berenices (Jeffries & James 99 data). Solid 
line is prediction of standard (corresponding to 
model Pleiades A) model of Pleiades composition. 
Small triangle on this line is IMq star of pleiades 
composition. Dashed-line is lithium fraction inter- 
polation in Pleiades. 
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Fig. 6. — Crosses : Li abundances for members 
of the Hyades (Thorburn 93 data). Triangles : 
Blanco I (Jeffries data). Lines show model com- 
putation for different masses at ^ 625 Myrs. Solid 
line is prediction of standard (corresponding to 
model Hyades A) model of Hyades composition. 
Dashed line is prediction of case E model (see sec- 
tion on opacities). Small diamonds on both lines 
represent IMq star. 
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said, predominant metals in opacity generation 
are iron, oxygen, neon, silicium, and magnesium. 
The full set of these metals represents more than 
70% of total opacity whereas any other accounted 
for metal do not contribute over 5%. In table 3 we 
give opacity contributions of 19 metals at BCZ for 
a IMq star extracted from monochromatic calcu- 
lations of Iglesias & Rogers (1996). Those results 
are obtained the following way : for every com- 
ponent wc compute total opacity for the stellar 
plasma mixing excluding that component. Differ- 
ence of opacity with the plasma including every 
components then allows to estimate to what level 
that element is responsible for opacity. 

Table 3: Repartition of elements in opacity cal- 
culations at BCZ for solar composition in typical 
main sequence (T ~ 2.10^Kp ~ 0.2g.cm~'^) and 
pre-main sequence (T ^ 4.10^Kp ~ 2g.cm~^) con- 
ditions. 



Element 


MS 


pre-MS 


H 


20% 


13% 


He 


11% 


6% 


C 


8% 


2% 


N 


5% 


1% 





■MM: 


19% 


Ne 


9% 


13% 


Mg 


3% 


11% 


Al 


< 1% 


1% 


Si 


3% 


12% 


S 


4% 


5% 


Ar 


1.5% 


<1% 


Ca 


2% 


< 1% 


Fe 


16% 


22% 


Ni 


1% 


2% 


For every following metal 






Na,P,Cl,K, 






Ti,Cr and Mn 


< 1% 


<1% 



We now try to determine the detailed composi- 
tion of the Hyades for relevant metals. Edvardsson 
et al.(1993) find a tight correlation between [0/Fe] 
and [Fe/H] in galactic disk stars, 
[O/Fe] = (-0.36±0.02)x[i^e/iJ]-(0.044±0.010) 



This relation once again suggests our Sun could 
be overabundant in oxygen when compared to 
other stars. Moreover Garcia Lopez et al. (1993) 
claim oxygen to hydrogen metallicity is close to 
[O/H] = -0.07 ± 0.05 dex in solar-like Hyades 
stars. This result is qualitatively consistent with 
the previous relation. If we apply this law, Hyades 
appear now metal-poor in oxygen compared to 
the Sun, this is in agreement with helium deter- 
minations from Ferryman et al. (1998). Neon 
is not detectable in solar-like stars photosphere 
but its ratio to oxygen measured in HIT regions 
is remarkably constant for different [O /H] (Meyer 
1989). Indeed neon, silicium, and magnesium are 
expected to vary like oxygen because of their com- 
mon origin in SNII. Consequently we decide to 
take the same variation for silicium and magne- 
sium than for oxygen. Let us just signal one caveat 
regarding silicium : Cayrel, Cayrel de Strobel, 
feCampbell 1985 find [Si/H] = 0.16 ± 0.05. Car- 
bon cannot be ignored at least on MS but its 
origins are still a matter of debate (Gustafsson et 
al. 1999). We let it vary like oxygen. 

We have generated opacity tables for this 
non-solar metal distributions using the devoted 
Lawrence Livermore National Laboratory website. 
By doing so we have chosen to model two cases. In 
the first one (case E), every metal varies like oxy- 
gen except iron-peak elements; in the second case, 
we let every metal vary like iron except the oxy- 
gen, to test the dependence of lithium depletion 
on oxygen (case F). In case E, which is the most 
realistic, lithium abundance is 2.6 dex (instead 0.9 
dex) for a IMq star. The shift of oxygen abun- 
dance from 0.127 to -0.07 dex increases lithium 
fraction by a large amount of 1.1 dex. Then the 
shifting of every metal but iron-peak ones to the 
level of oxygen increases again lithium fraction 
by 0.6 dex. This dramatic increase in lithium 
fraction is shown on figures 6 and 7. Lithium in 
solar-like Hyades stars becomes quite compatible 
with observations. There are two reasons for this. 
First, the decrease in metal fraction reduces BCZ 
temperatures and secondly it increases effective 
temperature for any given mass. Solar mass stars 
are shifted from ~ 5450 to ~ 5800 Kelvins. The 
large effect of metal might at first seem worry- 
ing. However one as to keep in mind two points. 
Firstly metals determine roughly 80 % of opacities 
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at BCZ during prc-MS. Secondly lithium depletion 
is a very steep function of temperature (see fig- 
ure 3). Furthermore we remark that such a large 
impact of mctallicity on depletion is reported in 
other recent studies (Chaboyer, Demarque, & Pin- 
sonneault 1995). Ventura et al. 1998 use OPAL 
opacity table although a slightly older version. 
They claim a decrease of 15% in the value of Z 
(around Z=0.02) makes the '^Li abundance to vary 
by almost two orders of magnitude. The present 
effect of oxygen (0.2 dex represents ~ 60%) which 
is half of metal fraction is quite compatible with 
this last result. 

We note that even lower mass stars appear in 
good agreement with observations. Figure 7 out- 
lines the importance of metals in lithium deple- 
tion. It also shows that oxygen plays a more de- 
terminant role than iron for pre-mainsequence evo- 
lution. 
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Fig. 7. — Predictions compared to observations 
for solar-like stars of Hyadcs composition. Solid 
line corresponds to Hyades' composition with ev- 
ery metal scaled on oxygen [O/H]=-0.07 dex ex- 
cept iron-peak elements (case E) . Dashed line cor- 
responds to model where all metal are scaled on 
iron [Fe/H]=0.127 dex (case A). For dashed dot 
line all metal arc scaled on iron except oxygon 
(case F) . Error-bars photospheric observed lithium 
are deduced from Thorburn et al. (1993) data 
for solar-type stars considering Te// within 100 K 
around 5450 K (case A) and around 5800 K (case 
E). 

If opacities are crucial at the base of convec- 



tion zone they are also of first importance in the 
atmosphere, as they determine depth where ma- 
terial becomes convective and this affects signifi- 
cantly the convection zone extension. The deeper 
convection starts in atmosphere the more effi- 
cient it is and the deeper it goes in the enve- 
lope. For low effective temperatures encountered 
on the Hayashi track many contributors to opac- 
ity have to be considered. The lowest tempera- 
ture we reach in computation is 3200 K and below 
5000 K molecules are not negligible anymore. It 
is therefore likely that improvements of present 
models may be made when using new atmosphere 
models (Hauschildt, AUard, & Baron 1999). Here 
we have mainly use Alexander & Ferguson opac- 
ity tables (Alexander & Ferguson 94) and checked 
that there was no important discrepancy induced 
when using low-temperatures Kurucz opacity ta- 
bles (Kurucz 1992). For case A composition star, 
final (and maximal) difference between otherwise 
similar models is less than 0.1 dex which is very 
small if we recall global depletion factor is more 
than 2 dex in this case. 

In typical pre-MS conditions H~ remains the 
main opacity source. The metals provide a larger 
part of free electrons when temperature lowers. 
Alexander & Ferguson opacity table have been 
computed for solar composition and is not adapted 
for case E. We can evaluate average ionisation 
states of all components in atmospheric thermody- 
namic conditions. For the coldest (~ 31647^") and 
most diffuse (1.5 10~^°g.cm~^) conditions, iron- 
peak elements and oxygen are respectively respon- 
sible for ~ 10 and 0.25 % of electron density. 
When opacity evaluation changes from Alexan- 
der to OPAL table, temperature is logT=3.75 and 
density is typically a few 10^'^ g.cm^'"*. We esti- 
mate iron-peak elements to be responsible for ~ 30 
% of electron density and contribution of oxygen 
being negligible. 

In case E metal fraction is reasonably scaled 
by oxygen which is the main mass contributor 
for heavy elements. Case E is therefore slightly 
warmer than case A model and its BCZ is at a 
lower temperature. At ~ 6 Myrs case A BCZ is 
located at O.3M0 where T = 3.99 10*^ if case E 
BCZ is located at O.4M0 where T = 3.72 lO^K. 
On the other hand by scaling every metal on oxy- 
gen we underestimate iron-peak element fraction 
(and effects) in the atmosphere where they could 
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provide up to 30 % of free electrons. If the right 
iron-peak element fraction was considered in the 
atmosphere, electron density would be enhanced. 
Convection would start closer to the surface and 
would therefore be less efficient. So we can expect 
that our values of BCZ temperature are in case 
E a bit too high. Figure 8 shows however that 
the changes among metal fractions are mainly felt 
at BCZ. Difference in temperatures as a function 
of mass never exceed 5% and changed opacities 
mainly affect BCZ but not the global structure. 



Age : 6 Myrs 
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Fig. 8. Relative temperature difference at a 
given depth between case E and case A models. 
Diamond and cross respectively locate BCZ of case 
A and E models. 

4. MACROSCOPIC EFFECTS: CON- 
VECTION AND ACCRETION 

4.1. Efficiency of convection 

The usual mixing-length theory parameter 
value relies on solar calibration. There is no reason 
to believe it is universal. Recent work based on 
hydrodynamical simulations (Ludwig et al. 1999) 
has investigated possible calibrations of mixing- 
length for solar- type stars. This work shows a 
dependence of a on effective temperature and sur- 
face gravity. For temperatures between 7100 and 
4300 K and surface gravities 2.54 <log g< 4.74 a 
varies from 1.3 to more than 1.7. Using this work 
we find negligible evolution of a on the Hayashi 
track for a 1 Mq object where surface gravity 
varies significantly but effective temperature re- 
mains almost constant. On the contrary once the 



radiative core becomes important and the star 
leaves the Hayashi track, the increase of the ef- 
fective temperature could justify a change of the 
a value. For a solar composition model this oc- 
curs at about 12 Myrs ie slightly before the end 
of the "^Li burning phase. In order to follow the 
effective temperature impact on a we have com- 
puted models of 1 Mq stars for Pleiades, solar, and 
Hyades composition. Effective temperature being 
sensitive to composition age and a parameters 
vary with both of them. Temperature differences 
on the Hayashi track are approximately 70 K if 
one changes from solar to Hyades composition. 
In region of interest in temperature and gravity, 
differences in a are not significant over such a 
narrow temperature range and the same a can be 
used. After leaving the Hayashi track tempera- 
ture differences increases up to 200 K but the star 
evolves towards solar main sequence conditions in 
a zone where there is a 'plateau' in a (Ludwig 
et al. 1999). Here again one can adopt a unique 
value so that for present study the change in o 
corresponds to temperature and not composition. 
On this point Ludwig et al.(1999) agrees with Fer- 
nandes et al. (1998) results who foimd that a is 
almost constant in the Sun and 4 low-mass stars 
systems having metallicities between solar and - 
0.31 dex, and helium content from 0.25 to 0.28. 
For these reasons we use the same values in the 
mixing length parameter whatever the composi- 
tion. Pleiades and solar composition being close, 
we only made distinction between this group and 
Hyades. Between these two groups composition 
determines the time spent on the Hayashi track. 
We have considered three different a values: on 
the Hayashi track, when the radiative stratifica- 
tion begins to influence surface temperature and 
on the main sequence where the value results from 
usual solar calibration. Following Ludwig et al. 
advices we have adopted a values they proposed 
after reestimate by 0.1 and use their values as 
scaling factors. If the star has solar or Pleiades 
composition, q:=1.935 before 15 Myrs ; a=1.85 
between 15 and 22 Myrs and a=1.766 after 22 
Myrs which is our solar-calibration main sequence 
value. If the star has Hyades composition, we 
use the same values but at different times and 
change the limits to 20 Myrs and 28 Myrs. Higher 
mctallicity results in slower contraction and evo- 
lution. In Hyades case only the initial value of a 
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is important because after 20 Myrs '^Li burning is 
over. 

Results of the hydrodynamical calibration must 
be considered cautiously. Firstly, because low 
temperatures molecular opacities are not included 
in Ludwig et al. calculations and we are presently 
exploring low effective temperatures regions. Sec- 
ondly a recent result about t Pegasi binary systems 
calibration predicts opposite evolution of a with 
effective temperature (Morel et al 2000). Thirdly 
2D calculations have to be confirmed by 3D hy- 
drodynamical calculations. 

In these calculations, ''iz depletion is increased 
(table 2, case B) as convection zone extends deeper 
to higher temperatures with more efficient convec- 
tion. We note nevertheless that such modification 
leads to less lithium destruction than using a full 
spectrum of turbulence as Ventura et al. (1998). 

4.2. Mass accretion 

Observations of T Tauri stars accretion lumi- 
nosities lead to mass accretion rates spanning from 
a few 10^*^ up to a few 10~^Mo/yr (Hartigan 
1995, Gullbring et al. 1998). This wide range 
probably originate from real star to star differ- 
ences altough such low accretion rates are very 
difficult to evaluate (Hartmann 1997) and suffer 
from large uncertainties. Accretion can affect "^Li 
stellar photosphere abundances in mainly three 
different ways. First, it has a structural impact 
as it modifies the stellar mass, consequently the 
gravitational potential then stratification change. 
Secondly, providing ISM abundance material to 
the surface of the star, has a direct chemical im- 
pact. Thirdly, accretion should also change stellar 
boundary conditions which is probably the most 
difficult part of the accretion phenomenon to mod- 
elizc. 

Our hydrostatic calculation takes accretion into 
account in a crude fashion. The accreted mass 
modifies only the external layers of the star and 
does not directly affect global surface boundary 
conditions (pressure, temperature, luminosity). 
Mass is simply added on the outermost layer of the 
star. We limit our study to accretion rates below 
lO~''M0/yr and following Hartmann (1997), we 
consider global accreted mass of a few 10~^ Mq. 
As recent hints suggest that accretion could last 
longer than usually believed and perhaps deal with 
larger accreted masses (MuzeroUe et al. 2000), so 



we will then also investigate the effects of O.IMq 
mass accreted. 

4.3. A low global accretion mass 

We first evaluate "^Li photosphere variations 
through accretion by simply considering different 
nonaccreting stellar masses. For a 2% accretion 
mass, if the structural effect is dominant, final 
minimal ''Li photosphere fraction is the one of 0.98 
Mq non-accreting star. This minimal fraction is 
1.92, to compare with 2.13 for a 1 Mq for a so- 
lar composition. Then the maximal fraction can 
also be evaluated, if we consider that the chem- 
ical mixing is dominant. This value is obtained 
in considering a IMq star but artificially increase 
''Li photosphere fraction by diluting in external 
CZ the ''Li mass contained in O.O2M0 of ISM ma- 
terial. The ''Li maximal fraction is 2.61 for a 1 

Mq. 

A real accreting star has a lower mass all the 
way through accretion phases and depletes both 
its initial lithium and lithium it receives from ac- 
cretion. So, in the following we consider accretion 
impact on ''Li burning for solar composition star 
and 2% Mq total accreted mass at varying tem- 
poral rates. Starting with a 0.98 Mq stellar ob- 
ject we consider a 'fast accretion' rate model of 
lO~^M0/yr during 2 Myrs, we get a lithium con- 
tent of 2.12 dcx and then a 'slow accretion' rate 
model lO~*M0/yr and lO~^M0/yr during 1 and 
10 Myrs respectively which gives a lithium content 
of 2.08 dcx. This kind of simulation docs not af- 
fect depletion by more than 0.05 dex. The general 
trend of accretion is to lower ''Li fraction so we 
can conclude that structural effects are predomi- 
nant over chemical effects. The effect is neither 
sufficient to recover the agreement with observa- 
tional ''Li fraction nor to explain the spread in 
young clusters. 

Accretion will have different consequences on 
lithium burning if it could last long enough ie af- 
ter major ''Li pre-MS depletion phase. Near IR 
excess of very low accretion rates of lO~^M0/yr 
are currently not dectected (Hartigan, Edwards, 
& Ghandour 1995). The simple evaluation we 
make just above illustrates that a variation of 
only very few percent of stellar mass could have 
a nonnegligible chemical impact because after 10 
Myrs external convection zone is 10 % or lower of 
the stellar mass. Then, an accretion rate as low 
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as 10 ^Mo/yr could significantly change ZAMS 
lithium surface fraction. 

4.4. Larger accretion rates applied to the 
Pleiades composition 

90% or more of stellar final mass is accreted 
during short (less than one myr) class and class 
I stages who extend up to a few 10^ years before 
classical T Tauri phase (Andre, Ward-Thomson, & 
Barsony 1999 and references therein). We consider 
here both different final masses from 0.9 to 1.1 Mq 
and different globally accreted masses ie 2 10~^ or 
1O~^M0. The stars now have Pleiades composi- 
tion. As we have seen in the solar case, structural 
effect of accretion is maximum when accretion is 
fast, we therefore restrict our computation to 'fast' 
accretion process: 10~^ or 5 10~^MQ.yTS~^ during 
the first Myr and then 10"^ or 5 lO'^Mg.yrs"^ 
from 1 to 11 Myrs. In these conditions, accre- 
tion still increases Li depletion and confirms so- 
lar trend. In table 4 we give lithium abundance 
relative to hydrogen as function of mass and with 
or without accretion in case of Pleiades. The lower 
the mass the higher the accretion impact. Accre- 
tion could therefore explain initial lithium disper- 
sion that also seem to increase towards low effec- 
tive temperatures. 



Table 4: Lithium surface abundances at 70 Myrs 
for accreting and non-accreting Pleiades composi- 
tion models. 



Mass 


I.IM© 


IMq 


0.9Mo 


Effective 








temperature 


6050 


5720 


5350 


Observed mean 








lithium abundance 


2.92 


2.70 


2.57 


No accretion 








models 


2.95 


2.58 


1.66 


Global accreted 








mass = O.O2M0 


2.94 


2.56 


1.59 


Global accreted 








mass = O.IMq 


2.86 


2.39 


1.24 



For a comparison, we estimate observational 
lithium abundance dispersion in the data set of 
Soderblom et al. (1993) we divide into 200 K bins. 
For each bin we then compute mean abundance 
value and dispersion. Between 4200 and 6000 K, 



Table 5: Evolution in Mq/uv of required accre- 
tion rates to reach lithium surface abundances at 
Pleiades age. 

Age (myrs) Accretion IMq Accretion O.9M0 



2 


6.3 10-« 


7.910-* 


3 


6.3 10-8 


7.910-8 


4 


4.010"^ 


6.310-8 


5 


3.110-^ 


4.010-8 


6 


1.6 10-s 


2.010-8 


7 


10-8 


1.210-8 


8 


6.310-^ 


7.910-9 


9 


4.0 10-9 


4.010-9 


10 


2.510-9 


2.510-9 


15 


1.610-1" 


2.510-1" 


20 


510-12 


2.5 10-" 



we find that dispersion varies from 0.65 to ^ 0.1 
dex exhibiting the well-known general trend to de- 
crease when tempature increases. Around Te// = 
5720 K (1 Mq) it is ~ 0.12 while around Te// = 
5350 K (0.9 Mq) it is - 0.32. As can be seen in 
table 4 these results arc qualitatively and quanti- 
tatively comparable to predicted dispersion result- 
ing from accretion of 10 % of Mq. We conclude 
that accretion could explain early main sequence 
dispersion. Results suggest also that the more a 
star accretes the more it depletes Li because once 
again in these cases structural effects dominate 
chemical effects. "^Li is not refreshed at a suffi- 
cient level to compensate for the additive burning 
due to the lower mass. 

We now evaluate the minimal accretion rate 
necessary to counteract mass effect for stellar 
masses of 0.9 and 1 Mq (table 5). Every mil- 
lion years photosphere "^Li fraction disminishes 
by a given amount and knowing convection zone 
mass and ISM ''Li fraction we compute a 'nomi- 
nal' mass accretion rate necessary to exactly com- 
pensate losses. One could however argue that 1 
Mq Pleiades have not kept initial ^Li fractions. 
They indeed depleted from ISM value ~ 3.2 down 
to ~ 2.8. Taking this into accoimt we compute 
a new accretion rate in the following way ; if a 
fraction a of 'nominal' accretion is provided every 
Myr only 1 — a of the quantity of depleted ''Li 
with no accretion will finally be depleted. We find 
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that new accretion rates should be lowered from 
'nominar ones by 0.32 dex and we provide results 
in table 5. Results do not really change within 
considered mass range. When compared to ob- 
servations of present star forming regions (Calvet 
& Gullbring 1998; Muzerolle et al. 2000) these 
accretion rates seem slightly too high. The re- 
quired accretion rates exceed strongest reported 
accretions observations by approximately 0.5 dex 
which brings us close to the observational upper 
limits. 

It is clear that conclusion on the accretion 
rates rciquires improved detection of mass accre- 
tion estimates to rates as low as 10~^MQ/yr for 
a large number of very young clusters, especially 
for low mass stars and corresponding photospheric 
lithium content. Therefore preceding results are 
mainly illustrative rather than conclusive. More- 
over, the accretion process has to be approached 
hydro dynamically and maybe in considering peri- 
odic phenomena, to take into account the modified 
boundary conditions. 

5. MACROSCOPIC EFFECTS: THE IN- 
FLUENCE OF THE ROTATION 

Regarding solar-like stars, open clusters "^Li 
abundances suggest a general depletion over main 
sequence which is in contradiction with standard 
evolution codes results. We will examine in this 
section the role of the internal rotation firstly on 
the pre-mainsequence structure and secondly on 
the possible tachocline mixing process inside radi- 
ation zone. 

5.1. Rotation structural effects 

Within the last decade many theoretical works 
have been led in the field of angular momen- 
tum transport in stellar interiors. Such works 
encounter difficulties to explain rotation veloc- 
ity evolution. For instance models invoking hy- 
drodynamic angular momentum transport (Pin- 
sonneault et al. 1989) predict strong differen- 
tial rotation in solar interior which is contradic- 
tory with helioseismology results. Models with 
angular momentum transport induced by internal 
magnetic fields (Keppens, MacGregor, & Char- 
bonneau 1995) can fit observations with a core- 
enveloppe coupling time of the order of 10 Myrs. 
They however encounter problems in reproducing 



the slowest rotators in open clusters and early MS 
rotational evolution. 

Core rotation could still significantly differ from 
surface rotation. Rotation rate is presently ob- 
served to be constant in most of the solar radiation 
zone and varies with latitude in solar convection 
zone. However it could have been a varying quan- 
tity in the initial stellar radiative and/or convec- 
tive interior. We investigate here purely structural 
effects of varying internal rotation and examine 
the consequences on effective temperature, lumi- 
nosity, and lithium content. The models we con- 
sider have Hyades composition where the stars be- 
low IMq all have experienced an initial rotational 
history increase until ZAMS followed by a strong 
decrease (see section 5.2). On the other hand they 
are still at the beginning of ZAMS so there is prob- 
ably no long term MS effects on lithium. This 
cluster therefore seems to us quite appropriate to 
study early rotational effects. We compute three 
models in the case of very short-lived circumstellar 
disk of 0.5 Myrs. The first model (SBR), assumes 
solid body rotation throughout the star. The sec- 
ond (DR), assumes complete decoupling between 
radiative and convective zones and the last model 
(FR) assumes solid body rotation in convection 
zone and every stellar mass shell keeps its initial 
angular momentum in radiative stratification part. 
If the first model corresponds to zero coupling time 
between convection and radiation zones, this cou- 
pling time is infinite in the latter ones. 

A well known structural effect of rotation is the 
decrease of the effective temperature and to a less 
extent of stellar luminosity (Sills, Pinsonneault, & 
Terndrup 2000). For the SBR stehar model we 
compute equatorial velocity of 74 km/s on ZAMS 
(~ 57 Myrs). This gives rise to a temperature de- 
crease of only 23 K in comparison with the nonro- 
tating model. Under analog conditions, the poly- 
nomial formula of Sills et al. gives a discrepancy 
of 15 K. The same formiila at maximal speed age 
(110 km/s for ^ 30Myrs) leads to 50 K where 
we find 65 K in difference. Differences never ex- 
ceed 0.03 dex in luminosity and 2 % in radius at 
given time. The expected consequence of such a 
behavior over convection is a larger extension and 
higher BCZ temperature and density responsible 
of the higher depletion rate in "^Li. Temperature 
at the BCZ increases strongly when the effective 
temperature lowers. This mass-lowering is effec- 
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tivc when the star reaches its maximal rotation 
rate around 30 Myrs. Yet before ^ IMyrs the 
situation is reversed and rotation models surpris- 
ingly exhibit lower BCZ temperature and density. 
A maximum difference in temperature of 2 10'*K 
appears at 3 Myrs. It seems indeed that early 
rotational effects over BCZ are opposite to what 
they become afterwards. This trend changes the 
temperature sensitivity on Li depletion rate. 

Quantitative variation of lithium fraction be- 
tween slowly (Tdisfe=3 Myrs) and rapidly {Tdisk 
=0.5 Myrs) rotating models is very small : 0.1 
dex (Table 6). No difference exists between SBR 
and DR models because the rotation is the same 
before 20 Myrs. In fact, the rotation rate, any- 
where in the star, is independent of any coupling 
time at least until 20 Myrs and at this age ^Li 
pre-MS depletion is over. There are two reasons 
for that. Firstly radiation zone stratification is 
still hardly under-adiabatic until 30 Myrs as we 
already explained in section 2.1. Secondly stellar 
wind loss represents only 12% of initial momentum 
at 20 Myrs (vs ~ 99% at solar age). Contraction 
of the whole star therefore stays homologuous and 
coupling-time between zones has no impact. Un- 
less initial conditions are different from solid ro- 
tation for the initial fully convective body there 
are no requirements to bother about angular mo- 
mentum loss or exchanges during early pre-MS Li 
burning phases. 

In the preceding lines we have been concerned 
by structural effects as far as they result from a 
correction to local gravity. We do not note any 
significant impact of these effects on the present 
problem. Now there arc others means by which 
rotation might affect stellar structure. Convective 
movements arc sensitive to rotation through cori- 
olis force. This effect is known for long in earthly 
atmospheric global circulation notably in the in- 
tertropical region. Siess & Livio 1997 suggest that 
convective cells could be twisted by rotation al- 
lowing therefore a smaller amit to the convective 
zone. As long as precise hydrodynamical compu- 
tation are led on that phenomenon it is however 
difficult to say much more about this solution. 

5.2. The turbulence at the base of the con- 
vection zone 

We consider now that the rotation induces 
an hydrodynamical instability in the tachocline 



layer at the top of the radiation zone (Spiegel & 
Zahn 1992). Such an instability has been studied 
by BTCZ99 to interpret the helioseismic results 
(Kosovichev et al. 1997) and could partly be at 
the origin of lithium destruction during the main- 
sequence. The introduction of a time dependent 
turbulent term in the diffusive equation permits a 
better agreement on solar photospheric light ele- 
ments between observations and models. BTCZ99 
remark that in the present Sun, tachocline mix- 
ing is thin enough not to deplete '^Be. However 
tachocline mixing depends on rotation and differ- 
ential rotation which is poorly known in pre-MS 
stars. 

We reexamine this process introducing a more re- 
alistic rotation history in pre-MS. A IM© model 
of Hyades composition and including tachocline 
mixing will not deplete more than 0.03 dex in ^Be 
from formation until ZAMS and none afterwards. 
Beryllium depletion does not occur although rapid 
rotation and high BCZ temperatures. It is the 
reason why we do not discuss '^Be in this study. 

Skumanich law (Skumanich 1972) has been 
used to infer rotational time evolution in main- 
sequence (as in BTCZ99). But such a law cannot 
be extend to young pre-MS stars which initially 
rotate slowly and then experience strong accelera- 
tion toward ZAMS. We reevaluate here tachocline 
mixing in using a more realistic rotation law that 
should apply over both pre-MS and MS. So we 
adopt the approach of Bouvier, Forestini, & Al- 
lain (1997, hereafter BFA97) to calculate different 
sets of rotation evolution for 1 Mq stars under the 
three following assumptions : (i) stars rotate as 
solid bodies, (ii) stars are locked to a defined an- 
gular velocity as long as they interact with initial 
surromiding disk, (iii) magnetic wind braking acts 
all along pre-MS and MS and depending on rota- 
tion speed produces varying angular momentum 
loss rate: 

(1) 

at lio ivio 

(2) 

Figure 9 presents corresponding equatorial ve- 
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Table 6: Impact of internal rotation model on lithium fraction in dex for solar mass and Hyades composition 
star 



TdiskinMyrs 


No rotation 


SBR 


DR 


FR 


Lithium at 7 Myrs 


2.06 


+0.11 


+0.11 


-0.01 


Lithium at 20 Myrs 


0.96 


+0.07 


+0.08 


-0.07 


Lithium at 0.7 Gyrs 


0.91 


+0.07 


+0.07 


-0.05 



locity evolution for solar mass and composition. 
Wsat and K are adjusted to fit observations of sur- 
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0.1 1.0 10.0 100.0 1000.0 
Age (Myrs) 

Fig. 9. — Equatorial velovity vs age for stars of so- 
lar mass and composition. Prom fastest to slowest 
rotators Tdisk is 0.5,3,10 Myrs. 

face rotation rates. Surface magnetic field in solar 

type stars seems to increase with f2 up to ^ lOfio 
and then saturates (Saar 1996). Saturation value 
determines a transition in the efficiency of the 
wind to brake the star which has consequences 
over early MS rotation rates. Following BFA97 
we adjust ujsat=^^ ^o- K determines the rotation 
rate at solar age. Recent hclioscismic measure- 
ments (Corbard et al. 1997) give internal solar 
rotation rates. A latitudinal dependence of rota- 
tion frequency is observed in convection zone, from 
460 nHz at equator down to 370 nHz at 60 degrees 
latitude. Radiative zone then experiences rigid ro- 
tation at the 430 nHz level and at least as deep as 
0.4 solar radius. In present solid rotation models, 
we adjust the rotational evolution laws so that the 
frequency is 430 nHz at solar age. To reach this 
rotation rate at actual solar age we adjust K to 



3.25 IQ'^^ g .cm? .s in the expression of the braking 
law. This value is a little larger than the one used 
by BFA97 (2.7 lO^'^g.cm^.s) as we tend towards a 
slower rotation velocity than surface equatorial : 
2.710-<5rrf.s-i instead of 2.9 lO-^rds'^ 

We then consider three different Tdisk dura- 
tions for star/disk coupling time : 0.5, 3, and 10 
Myrs. 0.5 Myrs corresponds to a star that ceases 
disk locking evolution on its birthline or a few 0.1 
Myrs after depending on mass and composition. 
Three Myrs is the median disk lifetime as esti- 
mated by BFA97. Finally 10 Myrs corresponds to 
a persistent disk. At this age only 10 to 30 % of 
young stars still show the IR and mm radio emis- 
sion expected if an optically thick disk is present 
(Strom 1995). The star is kept on a velocity of 
9.1 lQ~^rd.s~^ until it uncouples from disk. Then 
owing to gravitational contraction it accelerates up 
to 2.70 10-Vd.s-i and 1.75 lO-Vds"! after ~ 30 
Myrs for Tdisk equal 10 and 0.5 Myrs respectively. 
Such rotation rates correspond to equatorial ve- 
locities that span from 20 km/s to 120 km/s. Af- 
terwards the stars decelerate rapidly towards the 
same low velocity whatever disk lifetime. At the 
Hyades age the stellar equatorial velocity is in the 
narrow range from 4 to 6 km/s. 

Metallicity has direct impact on rotation evo- 
lution. Wc evaluate rotation history for three sets 
of mctallicities : Pleiades, Sun. and Hyades. From 
the former to the later composition contraction 
time towards ZAMS increases. Peak velocities 
arc reached later and are slightly lowered mainly 
due to higher moment of inertia and radius. A 
IMq pleiad reaches maximal velocity at 29 Myrs 
whereas its Hyades conterpart reaches it at 33 
Myrs. Moreover radius being a bit larger the wind 
braking affects more surface rotation. For a typi- 
cal disk lifetime of 3 Myrs we estimate peak veloc- 
ity to be 120 km/s and 113 km/s for Pleiades and 
Hyades compositions. These values are reached 
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around 29 and 33 Myrs respectively. 

Purely structural effects of rotation have been 
included the same way as in the preceeding sec- 
tion. Rotational centrifugal acceleration effects 
are added to gravity. This modifies the equation 
of hydrostatic equilibrium and radiative gradient. 
At radius r and angular velocity oj, the present 
stellar code integrates the 'mean' effect of rota- 
tion by substracting 2rw^/3 to gravitation. Then 
rotation velocities are also taken into account in a 
macroscopic rotationally induced diffusion coeffi- 
cient Dt (equation 14 in BTCZ99) and tachocline 
thickness d (equation 11 in BTCZ99) that deter- 
mines tachocline mixing efficiency. These coeffi- 
cients are assumed to follow the scaling laws (sug- 
gested by BTCZ99): 

i^T a 17°-75±0.25 
^^^(1.3±0.1)/4 

5.2.1. The solar case 

Rotation induces a strongly increased deple- 
tion of photospheric ''Li. Structural modifications 
do not play any significative role on pre-main se- 
quence as well as on main sequence unless the 
stellar disk lifetime is very short. For a 3 Myrs 
disk lifetime we note no differences exceeding 0.1 
dex in ''Li between a model that includes rota- 
tion structural changes and a model that does 
not. Stellar structures and evolutions are simi- 
lar. On the contrary rotationally induced mixing 
dramatically changes ''Li history. The 3 Myrs disk 
lifetime star experiences an enhanced destruction 
of ''Li during radiative core developpement phase 
when compared to non rotation model. ''Li abun- 
dance is lowered from ~ 2.1 dex without rotation 
down to ~ 1.8 dex. This increases ''Li depletion 
during pre-main sequence in comparison with pre- 
vious calculations using only the Skumanich law 
(BTCZ99). But the final result is not much al- 
tered. ''Li depletion depends on disk lifetime Tdisk- 
Variation is important between Tdisfc=0.5 and 3 
Myrs but more moderate for longer disk lifetimes 
(table 7, figure 10, and 11). 

Typical scatter between long and short lived 
disk is 0.5 dex and tends to grow slightly on main 
sequence. This does not agree with open-cluster 
observations that shows no dispersion in ''Li abun- 
dances from Pleiades to Hyades as can be seen on 



previous plots. On the other hand spread reap- 
pears in much older clusters as M67 (Jones, Fis- 
cher, & Soderblom 1999). With a photospheric 
lithium of 1.1 dex, the Sun is expected to have ex- 
perienced a long-lived disk stars (figure 10), if this 
approach is correct. However one has to be very 
extremely caiitious when drawing such conclusion. 
Pre-MS depletion phase is overestimated as young 
open-clusters observations suggest. 

5.2.2. The Hyades case 

Tachocline thickness and macroscopic diffusion 
coefficient depends on the rotation history and 
the mctallicity. However if scaling laws of equa- 
tions 3 and 4 are correct, we compute that differ- 
ences in metallicity of 0.1 dex do not induce large 
enough modifications in rotation speed to signif- 
icantly change tachocline diffusion coefficient or 
thickness in pre-MS. The general behaviour no- 
ticed for the Sun is also observed for this cluster. 
But the computation shows that higher velocity 
rotating stars do exhibit larger ''Li depletion, this 
does not agree with present observations, indicat- 
ing that rapid rotators have higher ''Li rates as 
their slower counterparts (Soderblom et al. 1993 ; 
Garcia Lopez, Rebolo, & Martin 1994). This point 
suggests that tachocline mixing, as it is now intro- 
duced in calculations, and probably any mecha- 
nism that would mix radiatively stabilized layers 
on pre-MS, could be partly inhibited. It is note- 
worthy that low mass star rotation are presum- 
ably evolving slower. A 0.8 Mq star should reach 
peak velocity at 50 Myrs. Moreover such a star 
will retain its angular momentum longer so that 
tachocline mixing retains longer efficiency. This 
brings the star the wrong way for it is the lighter 
stars that present the largest discrepancy in ''Li 
content between current models and observations. 

6. SUMMARY AND PERSPECTIVES 

This study confirms that lithium offers an ex- 
tremely interesting insight over solar-like stellar 
structure and evolution due to its low depletion 
temperature. There is presently a discrepancy be- 
tween computations of classical stellar models and 
observations. Classical models predict no lithium 
depletion on main sequence whereas it seems to 
be observed in open-cluster lifetimes. It is rea- 
sonable to think that lithium evolution on main 
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Table 7: Lithium fraction in dex for solar mass and composition star as a frmction of circumstellar disk 
lifetime. 



TdiskinMyrs 


0.5 


3 


10 


20 


Lithium at ZAMS 


1.34 


1.80 


1.90 


1.92 


Lithium at 0.7 Gyrs 


0.84 


1.43 


1.56 


1.61 


Lithium at 4.6 Gyrs 


0.35 


0.95 


1.10 


1.16 
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Fig. 10. — "^Li surface abundances in solar mass 
and composition stars. From the lowest to the 
highest depletion rate : first model only includes 
microscopic diffusion and three following models 
add tachocline mixing with disk coupling time 
10, 3, and 0.5 Myrs respectively. Also plotted 
M67 and solar surface lithium abundances. Ev- 
ery model exhibits the same effective temperature 
of 5777 K at 4.6 Gyrs (table 2 : solar models). 

sequence is connected to some slow mixing rota- 
tionally induced proccss(cs) that occurs at the top 
of the radiative zone. In the specific case of our 
Sun such a process is supported by differences be- 
tween theoretical and measured sound speed in 
this part of the star (BTCZ 1999) and by the 
agreement between photospheric observations and 
predictions for helium, metals, and lithium. Dur- 
ing pre-mainsequence, the lithium problem is re- 
versed and classical uptodate theoretic;al models 
generally predict too strong depletion. We confirm 
that this depletion strongly depends on metallic- 
ity, in the calculation. It is not evident that it 
is currently observed in young open-clusters. In 
this regard Coma Ber and Blanco I rather suggest 
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Fig. 11. — "^Li surface abundances in solar mass 
and Hyades case E composition stars. Same pre- 
sentation than above. Also plotted Blanco I and 
Hyades surface lithium abundances. Microscopic 
diffusion model has an effective temperature of 
5802 K at 625 myrs. At this age all other models 
exhibit very similar effective temperatures (from 
5800 to 5802 K) 

an age dependence. Moreover, early dispersion in 
lithium abundances has to be explained. 

In this paper, we have studied the impact of 
different microscopic or macroscopic processes on 
pre-MS lithium depletion. First we note the im- 
portant role of the choice of meshes and time steps 
in the very early phase (before 10 Myrs). Sec- 
ondly, we remark that the interpretation of data is 
presently difficult as Li depiction in open-clusters 
is a strong function of effective temperature and 
it is therefore very difficult to reliably relate an 
observed lithium fraction to a given stellar mass 
just because this effective temperature depends on 
detailed metal fractions that still have to be accu- 
rately determined (see figure 6 compared to figures 
7 and 11). 
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Microscopic processes are mainly dominated by 
opacity coefficients and therefore through metal 
and helium fractions. We experience in models 
that the level of sensitivity to composition (espe- 
cially on oxygen and iron but not only) proves how 
difficult the present situation is. This suggests 
that the dispersion in lithium abundances could 
be related to small dispersion in metal fractions or 
among metal fractions. This is all the more true 
as part of the observed scatter is surely related to 
observational questions. Rotation history might 
induce such difference but on the other hand King, 
Krishnamuthi , & Pinsonneault 2000 note the "^Li 
abundance determination to be correlated to (and 
probably biased by) chromospheric activity. Gen- 
uine dispersion could also be related to larger he- 
lium fraction variations which can not be ruled out 
even over small scales in the ISM (Wilson & Rood 
94 and references therein). In the solar case, where 
all the abundances are known within 10%, we get 
the good order of magnitude of lithium depletion 
in adding a macroscopic effect to the implicit mi- 
croscopic effect in pre-mainsequence. We do not 
observe also a large problem in Hyades when using 
observed oxygen and iron are measured. Neverthe- 
less there is a contradiction between the large de- 
pendence of lithium depletion in theoretical mod- 
els in comparison to the apparent quasi metal- 
licity independence of the abundances in young 
open-clusters. To solve this problem we encour- 
age measurements of photosphcric abundances as 
carbon, oxygen, silicium, and iron on young T- 
Tauri stars. Similarly helium fraction needs to be 
carefully checked within clusters and if possible 
from star to star. Such additional data can be de- 
duced from binary systems, high mass stars, and 
importantly in solar-like stars thanks to the future 
asteroseismology developments. 

Macroscopic processes are at play through ac- 
cretion, rotation, and convection modelling. Such 
processes are now crucial to understand accurately 
stellar evolution (Pinsonneault 1997). Accretion 
is likely to explain dispersion in Li fraction mea- 
surements, but it fails to bring mean abundances 
to an acceptable level. In Pleiades case, required 
accretion rates at a given time arc at the upper 
limit of observations for young T Tauri. To im- 
prove significantly the situation they should last 
so long that total accreted mass would indeed be 
~ 20 — 30% of solar mass during the first Myrs fol- 



lowing star formation. On the other hand it is also 
possible for young stars not to accrete from classi- 
cal gas & dust disks but from planetesimals of pro- 
toplanetary systems which are not detectable to- 
day. Such material infall could supply lithium and 
others metals in low mass convective enveloppes of 
post-TTauri stars. 

Rotationally induced mixing gives satisfactory 
results regarding photospheric solar lithium abun- 
dance. These results are moreover consistent with 
Hyades open-cluster observation when oxygen and 
iron fractions are correctly ie separatedly consid- 
ered. Rotation mixing process we consider is how- 
ever unable to explain initial spread in "^Li surface 
abundances. Treating more properly the fruit- 
ful MS tachocline mixing (BTCZ 1999) in pre- 
MS leads to results not supported by observations 
(King, Krishnamurthi, & Pinsonneault 2000 and 
references therein) : we find a rotation-^Lz deple- 
tion correlation on ZAMS. However it is important 
to keep in mind two points. Firstly we recall that 
the present rotationally induced mixing is not re- 
lated to angular momentum transport in the stars. 
It therefore offers a possible but certainly partial 
vision of lithium MS evolution. ZAMS rotational 
velocity should affect ^Li depletion because of sub- 
sequent angular momentum redistribution (Pin- 
sonneault et al. 1999 and references therein). Sec- 
ondly the rotationally induced mixing process we 
invoke has necessary a transitory regime that we 
did not take into account. For these two reasons 
"^Li depletion predictions of the tachocline mixing 
are more reliable between the Hyades and older 
clusters than absolute predictions from star for- 
mation. Hyades solar-like stars are already slow 
rotators and might have lost most of their angular 
momenta. From Hyades to older clusters, lithium 
should traduce long term (stationary) mixing ef- 
fects, and indeed the tachocline mixing give rise 
to the observed characteristic time (~ IGyr) for 
^Li depletion on MS. 

Finally it is possible that the convection pa- 
rameter amit is sensitive to rotation and hydrody- 
namical simulations suggest it should be changed 
along Hayashi track. The effect we mention is not 
very large, contrary to some previous studies but 
more work is needed to be sure. We have em- 
phasized the absence of a well established radia- 
tive stratification in central parts during phases 
preceeding 20 Myrs which is a tricky point. The 
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first phase of lithium burning occurs at the fron- 
tier between cooler convective medium and hot- 
ter radiative (but almost convective) medium that 
lies beneath. This means that very slight addi- 
tional perturbation (through overshooting or in- 
stability for instance) could give rise to a much 
stronger lithium depletion. This also moans that 
very slight stabilisation phenomenon, related to 
magnetic field as it is suggested by Ventura et al. 
(1998) could save a large fraction of lithium. 

All these remarks on macroscopic phenomena 
encourage complementary studies implying hydro- 
dynamics approach. This is also true to check if 
the star is initially fully convective and rotates as 
a solid body. The answers to those questions can 
only be extracted from earlier phases that include 
protostellar collapse and take into account hydro- 
dynamical processes such as accretion in a consis- 
tent and more rigourous manner. 
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